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THE  SHOCK- INDUCED  FOIARIZATION  OF  DIELECTRICS 

ABSTRACT 

In  previous  experiments  it  was  shown  that  a  displacement  current  is 
generated  in  suitable  dielectric  materials  subjected  to  explosively  produced 
shock  waves.  A  theoretical  description  of  the  shock- induced  displacement 
current,  which  takes  into  account  the  electrical  properties  of  the 
dielectric  and  the  shock  propagation  characteristics  of  the  media,  is 
used  to  interpret  data  obtained  from  experiments  employing  polystyrene 
as  the  dielectric.  The  theory  predicts  that  the  initial  current  is 
directly  proportional  to  the  area,  inversely  proportional  to  the  thickness 
of  the  dielectric,  and  independent  of  the  resistance,  all  of  which  agree 
with  earlier  experimental  findings.  The  data  show  that  most  of  the  Increase 
in  the  dielectric  constant  is  due  to  an  increase  in  the  density  of  the 
specimen,  with  a  rise  in  temperature  playing  a  lesser  role. 
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INTRODUCTION 


In  recent  experiments,  Q.  E.  Hauver '  showed  that  a  displacement  current  is 
generated  in  a  suitable  dielectric  material  subjected  to  an  explosively  produced 
shock  wave.  The  existence  of  the  displacement  current  was  established  by 
placing  the  dielectric  between  conducting  electrodes,  one  of  which  was  in  direct 
contact  with  an  explosive  charge.  When  the  charge  was  detonated,  a  shock  was 
transmitted  into  the  dielectric  material  through  the  electrode  next  to  the  charge. 
The  peak  pressure  was  controlled  by  varying  the  thickness  of  the  intervening 
electrode  or  buffer  plate.  As  the  shock  entered  the  dielectric,  an  electric 
current  was  detected  in  a  resistor  connecting  the  two  electrodes.  Since  there 
were  no  other  sources  of  EMF  in  the  circuit,  the  current  in  the  resistor  was 
proportional  to  a  shock- induced  displacement  current  in  the  dielectric. 

It  has  been  postulated^"  that  the  shock  wave  produces  a  net  polarization  of 

the  dielectric  by  partially  orienting  the  elementary  dipoles  of  which  the  material 

is  composed.  It  is  now  believed  that  preferential  elongation  or  compression  of 

the  dipoles  depending  upon  their  orientation  is  a  more  plausible  physical 

mechanism.  A  difference  in  the  mass  associated  with  the  ends  of  the  polar 
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molecules  has  been  suggested  as  the  most  likely  reason  for  the  shock  front  to 
produce  a  preferential  displacement  of  the  positive  and  negative  charges. 

The  research  to  be  discussed  in  this  report  was  conducted  for  the  purpose 
of  establishing  the  relation  between  the  electric  current  generated,  the  electrlca 
properties  of  the  dielectric,  and  the  shock  propagation  characteristics  of  the 
media.  A  theoretical  description  of  the  shock- induced  displacement  current  is 
derived  and  used  to  interpret  data  obtained  from  experiments  employing 
polystyrene  as  the  dielectric  material.  As  yet,  results  of  the  work  have  not 
been  interpreted  in  terms  of  molecular  mechanics  describing  the  interaction  of 
the  shock  front  with  the  polar  molecules  of  which  the  dielectric  is  composed. 

THEORETICAL  DESCRIPTION 

A  mathematical  description  of  the  shock-induced  displacement  current  is 
based  on  the  physical  model  illustrated  in  Figure  1(a).  At  a  time  t,  the  shock 
will  have  moved  into  the  dielectric  a  distance  Ut,  and  the  rear  surface  will 
have  moved  forward  a  distance  ut.  Ahead  of  the  shock,  the  electric  properties 
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are  described  by  a  dielectric  constant  K*  an  electric  displacement  D,  an 
electric  polarization  P,  and  an  electric  field  strength  E.  It  is  assumed 
that  the  dielectric  is  isotropic  (D  ■  +  P)  and  that  all  vectors  are  oriented 

perpendicular  to  the  shock  front.  Primed  quantities  are  used  to  describe  the 
electric  properties  behind  the  shock  front,  it  being  assumed  that  the  properties 
are  constant  throughout  the  shocked  region.  Across  any  boundary  between  two 
dielectric  media,  the  normal  component  of  the  electric  displacement  is 
continuous.  Since  the  electric  displacement  vector  is  perpendicular  to  the 
shock  front 

D-D'.  (1) 

Within  a  parallel-plate  condenser, 

D  -  Q/A,  (2) 

where  Q  is  the  charge  on  the  condenser  and  A  is  the  area  of  the  plates. 

Ahead  of  the  shock,  the  electric  displacement  is  related  to  the  electric 
field  strength  by  the  relation 

D  -  K«0E,  (3) 

where  cQ  *  8.85  x  10"  farad/m.  Behind  the  shock  front,  it  is  postulated 

that  the  polarization  consists  of  an  electrical  contribution  P*  and  a 

mechanical  contribution  Ps  resulting  from  the  shock  wave.  The  electrical 

contribution  to  the  polarization  is  proportional  to  the  local  electric  field 

strength;  i.e.,  P'  -  €0X'E'  where  X'  is  the  para-electric  susceptibility.  On 

the  basis  of  the  preceding  assumptions  the  electric  displacement  behind  the 

shock  is  given  by  D'  -  c  E'  +  P'  +  P  or 

o  s 

D*  -  K'c0B«  +  Ps,  .  (4) 

where  K'  *=  (l  +  X* )  and  P_  is  the  net  dipole  moment  per  unit  volume  produced 

B 

by  the  shock  wave.  When  equations  (1)  through  (4)  are  solved  for  the  electric 
field  strength,  one  obtains 

B  -  Q/K«0A, 

B»  -  Q/K'«0A  -  P0/K'€o 
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(5) 

(6) 


An  examination  of  Figure  l(a)  shows  that  the  potential  drop  between  the  two 
plates  la  given  by 

V(t)  =  E(Xq  -  Ut)  +  E‘  (U  -  u)t.  (7> 

Since  Kirehoff’s  law  requires  that  V(t)  +  IR  =  0,  one  obtains  the  following 
equation  for  Q: 

IR  +  (Q/K€0A)  (Xq  -  Ut)  +  (Q/K'eoA  -  Ps/K'eo)  (U  -  u)t  =  0,  (8) 

after  substituting  the  quantities  for  E  and  E'  given  by  equations  (5)  and  (6). 


It  is  convenient  to  represent  the  behavior  of  the  dielectric  by  the 
equivalent  circuit  shown  in  Figure  l(b).  The  current  in  the  equivalent  circuit 
is  governed  by  the  equation 

Q/C1  +  Q/C2  +  IR  «  £.  (9) 

An  examination  of  equations  (8)  and  (9)  shows  that 


|  =  Ps  (U  -  u)t/K«eo,  (ID) 

1/C1  =  (U  -  u)t/K*eoA,  (11) 

i/c2  =  (XQ  -  Ut)/KeoA.  (12) 

It  is  noted  that  the  electromotive  force  is  a  saw-tooth  wave  having  a  period 
equal  to  the  transit  time  of  the  shock  through  the  dielectric.  If  the  time 
constant  of  the  circuit  is  much  less  than  the  period  of  the  electromotive  force; 
i.e. ,  RC  «  XQ/U,  it  may  be  assumed  that  the  entire  potential  drop  appears 
across  the  capacitance.  In  this  approximation  equation  (8)  reduces  to 


(Q/KeoA)  (XQ  -  Ut)  +  (<4/K'€0A  -  Ps/K»eo)  (U  -  u)t  -  0.  (15) 

The  Initial  current  can  be  obtained  by  differentiating  equation  (13)  and 
evaluating  the  derivative  for  t  and  Q  equal  to  zero.  It  is  found  that 

IQ  -  (dft/dt)t  b  Q  -  Ps  (U  -  u)KA/K'X0.  (14) 

Equation  ('14)  predicts  that  the  initial  current  is  directly  proportional  to  the 
area,  inversely  proportional  to  the  thickness  of  the  dielectric,  and  ie 

independent  of  the  resistance,  all  of  which  agree  with  earlier  experimental 
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findings  .  The  experiments  to  be  described  in  this  report  were  conducted  in  an 
effort  to  confirm  the  theoretical  model  and  evaluate  the  constants  K’  and  Pg  for 
polystyrene. 
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KXfBRDflBBAL  APBftRAIUS 


The  experimental  apparatus  used  to  study  the  shock- Induced  polarisation 
in  polystyrene  is  shown  schematically  in  Figure  2.  An  explosive  charge  is 
used  to  produce  a  plane  shock  front  in  an  aluminum  buffer  plate  located 
between  the  explosive  and  the  dielectric  specimen.  The  thickness  of  the 
buffer  plate  is  used  to  control  the  peak  pressure  to  which  the  specimen  is 
subjected.  In  order  to  remove  electrostatic  effects  associated  with  the 
ionized  detonation  products,  the  buffer  plate  is  maintained  at  ground  potential. 
The  dielectric  is  Inserted  between  the  buffer  plate  and  a  second  electrode  which 
has  the  same  diameter  as  the  specimen.  The  specimen  and  second  electrode  are 
enclosed  in  a  Teflon  receptacle,  which  is  filled  with  carbon  tetrachloride.  The 
carbon  tetrachloride  reduces  the  impedance  mismatch  at  the  lateral  surfaces  of 
the  dielectric  thereby  reducing  the  rarefaction  wave  that  would  otherwise  result 
from  a  free  boundary.  The  buffer  plate  is  connected  to  the  second  electrode 
through  a  resistor.  The  shock-induced  displacement  current  is  observed  by 
monitoring  the  voltage  across  the  resistor  using  the  upper  beam  of  a  type  555 
Tektronix  dual-beam  oscilloscope.  A  representative  voltage-time  signal  is 
shown  in  Figure  3.  The  duration  of  the  signal  corresponds  to  the  transit  tins 
of  the  shock  through  the  specimen.  If  the  specimen  is  sufficiently  thin,  the 
attenuation  of  the  wave  while  traversing  the  specimen  will  be  negligible,  and 
one  can  assume  the  shock  velocity  to  be  constant. 

The  shock  parameters  are  determined  experimentally  by  measuring  the  free- 
surface  velocity  of  the  aluminum  buffer  plate  near  the  dielectric  specimen. 

The  arrival  times  of  the  free  surface  at  four  locations  in  space  are  detected 
by  a  set  of  probes  connected  to  a  suitable  pulse-forming  network  and  the  lower 
beam  of  the  dual-beam  oscilloscope.  The  arrival  times  are  indicated  by  the 
rapid  vertical  deflections  of  the  lower  trace  illustrated  in  Figure  3* 

DST2RMIHATI0H  OF  SHOCK  HUESSURB 

The  pressure  produced  in  the  dielectric  material  is  determined  by  a 
method  similar  to  that  outlined  by  J.  X.  Walsh  and  co-workers^,  The  particle 
velocity,  u,  resulting  from  the  shock  in  the  buffer  plate  is  computed  from  the 
approximate  relation 

0-5) 
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where  u.  is  the  free-surfaee  velocity.  The  shock  pressure  in  the  buffer 

X  S 

plate  is  determined  from  the  p  -  u  representation  for  the  Hugoniot  data  of 
2024  T3  aluminum.  As  shovn  in  Figure  4,  the  curve  for  aluminum  is  reflected 
about  a  vertical  line  through  the  point  representing  the  pressure  and  particle 
velocity  existing  in  the  buffer  plate.  The  intersection  of  the  p  -  u  curve 
representing  the  Hugoniot  data  for  polystyrene  with  the  reflected  p  -  u  curve 
for  aluminum  defines  the  pressure  and  particle  velocity  existing  behind  the 
shock  transmitted  across  the  interface  into  the  polystyrene.  The  shock  velocity 
in  the  polystyrene  is  now  determined  from  the  conservation  of  momentum: 

P  =  pQUu.  (16) 

The  shock  velocity  determined  from  equation  (16)  for  the  l/8  in.  polystyrene 

■z 

apecimen  was  5-26  x  Hr  m/ sec .  The  shock  velocity  determined  from  Figure  3 
by  taking  the  ratio  of  the  thickness  to  the  transit  time  measured  from  the 

z 

pulse  duration  is  also  5*26  x  l(r m/sec. 

DETERMINATION  OF  CIRCUIT  CAPACITANCE 

In  deriving  equation  (13),  it  was  assumed  that  RC  «  XQ/U.  To  check  the 
validity  of  this  assumption,  and  consequently  the  validity  of  equation  (13), 
it  is  necessary  to  measure  the  capacitance  of  the  circuit.  When  an  inductance, 

L,  is  inserted  in  series  with  the  resistor,  the  current  oscillates  with  a 
frequency  whose  period,  T,  is  given  by 

T  =  2jt  VEc.  (17) 

A  typical  trace  obtained  from  an  experiment  using  an  inductance  is  shown  in 

Figure  5(a)*  The  time  between  the  first  two  peaks  was  measured  to  determine 

the  period  associated  with  each  value  of  the  inductance.  As  shown  in  Figure  5(b), 

the  capacitance  can  be  obtained  by  plotting  the  linear  relation  between  T^  and  L, 

2 

the  slope  of  the  line  being  4n  C. 

*12 

The  data  shown  in  Figure  5(b)  yields  a  capacitance  of  8.8  x  10  farads 
for  a  1/8  in.  thick  polystyrene  specimen  having  a  1.0  in.  diameter.  The 
associated  RC  time  constant  is  8.2  x  lcf 10  sec.  for  a  resistance  of  95  ohms. 

On  the  other  hand,  XQ/U  is  6.0  x  10“^  sec.,  which  readily  satisfies  the 
condition  RC  «  XQ/U. 
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DETERMINATION  OF  DIELECTRIC  FROFERTEBS 

The  displacement  current  was  carefully  measured  for  a  series  of  four  tests 
using  1  In.  diameter  polystyrene  specimens.  The  experimental  conditions 
associated  with  each  test  are  sumnarlzed  below: 


Test 

Specimen 

Thickness 

Resistance 

Pressure 
(from  ufB) 

1 

1/8  In. 

32  ohms 

108  kllobars 

2 

1/8  in. 

93  ohms 

110  kllobars 

3 

3/l6  in. 

93  ohms 

110  kllobars 

4 

1/4  In. 

93  ohms 

110  kllobars 

As  long  as  the  resistance  Is  sufficiently  small,  compared  with  the  reactance  of 
the  stray  capacitance,  the  charge  on  the  stray  capacitance  can  be  neglected; 
i.e.,  dQ/dt  »  I  where  Q  Is  the  charge  on  plates  surrounding  the  specimen.  For 
each  test  the  current-time  curves  were  integrated  to  obtain  Q  as  a  function  of 
time;  the  resulting  curves  are  shown  In  Figure  6. 

To  facilitate  the  determination  of  K’  and  Pq  equation  (13)  was  put  into  the 
following  form: 

F(Q)  -  Q(X0  -  Ut)/KeoA(U  -  u)t  -  P8/K'eo  -  Q/K»eoA.  (18) 

It  is  noted  that  F(Q)  Is  a  linear  function  of  Q  whose  slope  is-l/K'€0A  and 
whose  Intercept  is  P  /K’e  .  Values  for  F(Q)  were  computed  from  the  Q(t)  curves 
presented  In  Figure  6  using  experimentally  determined  values  for  U  and  u.  Curves 
representing  the  function  F(Q)  are  presented  In  Figure  7-  It  is  noted  that  a 
linear  relation  exists  between  F(Q)  and  Q  for  both  tests  using  l/8  In.  thick 
polystyrene  specimens.  An  Increasing  departure  from  linearity  exists  for  the 
3/16  In.  and  1/4  In.  specimens  respectively.  For  the  thicker  specimens.  It  is 
possible  that  the  low-pressure  tall  following  the  shock  results  in  substantial 
variations  In  the  dielectric  properties  throughout  the  shocked  region,  which  is 
contrary  to  the  assumptions  used  In  deriving  equation  (l B).  For  the  l/8  In. 
specimens,  the  data  follow  the  linear  relation  between  F(Q)  and  Q  predicted  by 
the  theory,  and  accurate  values  of  K*  and  Ps  can  be  obtained  from  the  slopes  and 
Intercepts  respectively.  The  results  are  summarized  on  the  next  page. 
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Test 


p 

s  2 
(coul/m  ) 

K* 

X* 

X'/X 

p'/p 

3.59  x  10'5 

3.29 

2.29 

1.5 

1.6 

3.75  x  10"5 

3.29 

2.29 

1.5 

1.6 

The  susceptibility  of  the  polystyrene  behind  the  shock  was  computed  from  the 
dielectric  constant  using  the  relation  X*  =  K*  -  1.  For  the  unshocked 
material  the  dielectric  constant  is  2.55;  and  the  corresponding  susceptibility, 
1.55.  Since  shock  compression  of  the  polystyrene  crowds  more  dipoles  into 
a  unit  volume,  the  ratio  X*/X  should  be  approximately  equal  to  p'/p  provided 
the  shock  has  not  greatly  increased  the  temperature  of  the  specimen.  An 
increase  in  temperature  reduces  the  para-electric  susceptibility  and 
consequently  the  ratio  X*/X  is  expected  to  be  less  than  p'/p.  However,  it  is 
interesting  to  note  that  the  ratio  X'/X  is  only  Ojo  less  than  the  ratio  p*/p 
indicating  that  most  of  the  change  in  susceptibility  is  accounted  for  by  simple 
compression  of  the  specimen. 


SUMMARY 

A  theory  has  been  developed  which  relates  the  shock-induced  displacement 

current  to  the  electrical  properties  and  shock  propagation  characteristics  of 

the  dielectric.  The  theory  predicts  that  the  initial  current  is  directly 

proportional  to  the  area,  inversely  proportional  to  the  thickness  of  the 

dielectric,  and  independent  of  the  resistance,  all  of  which  agree  with  earlier 

experimental  findings.  The  response  of  the  dielectric  to  a  shock  wave  is 

described  in  terms  of  a  function  F(Q),  which,  according  to  the  theory,  is  a 

linear  function  of  Q.  Empirical  curves  representing  F(Q)  were  obtained  for 

polystyrene  specimens.  For  sufficiently  thin  specimens,  F(Q)  was  found  to  be 

a  linear  function  of  Q  and  the  theory  was  used  to  evaluate  P  and  K'  for 

s 

polystyrene  subjected  to  a  110  kilobar  shock.  The  data  show  that  most  of  the 
increase  in  the  dielectric  constant  produced  by  a  shock  wave  is  due  to  the 
increased  density  of  the  specimen  with  an  increase  in  temperature  playing  a 
lesser  role. 
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Fit.  I- SIMIAN  ILLIISTIATMt  (a)  THE  PHYSICAL  NOBEL  NUB  TO  DEKIIIE 
THE  NET  POLARIZATION  PRODUCED  BY  A  SHOCK  KAYE  IN  A  BIELECTIIC, 
All  (bl  YNE  EQUIVALENT  CiCNIT  lESBLTM  PICK  AN  ANALYSIS 
BF  TNE  NOBEL.  C|  REPREKITS  TIE  STOAT  CAPACITANCE  AISBCIATO 
VITN  THE  CWCBIT. 


FIGURE  2:  EXPERIMENTAL  APPARATUS  USED  TO 
STUDY  THE  SHOCK- INDUCED  POLARIZATION  IN 
P0LY8TYRENE. 
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FIGURE  3. -TYPICAL  OSCILLOGRAM  OBTAINED  FROM  A  TEST  USING 
A  1/8  INCH  POLYSTYRENE  SPECIMEN.  SINCE  THE  SPECIMEN 
ACTS  AS  A  CONSTANT  CURRENT  GENERATOR,  THE  VERTICAL 
DISPLACEMENT  IS  PROPORTIONAL  TO  THE  LOAD  RESISTOR.  THE 
ARRIVAL  TIMES  OF  THE  FREE  SURFACE  AT  FOUR  LOCATIONS  IN 
SPACE  ARE  DETECTED  BY  A  SET  OF  PROBES  CONNECTED  TO  A 
SUITABLE  PULSE- FORMING  NETWORK  AND  THE  LOWER  BEAM  OF 
A  DUAL- BEAM  OSCILLOSCOPE.  THE  ARRIVAL  TIMES  ARE  INDICAT¬ 
ED  BY  THE  RAPID  VERTICAL  DEFLECTIONS  OF  THE  LOWER  TRACE. 
THE  FREQUENCY  OF  THE  REFERENCE  SINE  WAVE  IS  50  MC  FOR 
THE  UPPER  BEAM  AND  10  MC  FOR  THE  LOWER  BEAM.  THE 
TRACES  ARE  DISPLACED  ALONG  THE  TIME  AXIS  BECAUSE  THE 
LOWER  BEAM  IS  TRIGGERED  ABOUT  0.2/iSEC  LATER  THAN  THE 
UPPER  BEAM. 
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TRANSMITTED  INTO  THE  POLYSTYRENE 


FIGURE  A 


FIGURE  B 

FIGURE  5 -(a)  TYPICAL  OSCILLOGRAM  OBTAINED  BY  PLACING  AN 
INDUCTANCE  L  IN  SERIES  WITH  THE  LOAD  RESISTOR.  THE  INIT¬ 
IAL  CAPACITANCE  IN  THE  CIRCUIT  IS  OBTAINED  FROM  THE  PERIOD 
OF  THE  FIRST  OSCILLATION.  SINCE  THE  PERIOD  IS  GIVEN  BY 
T*  2tt^/LC,  T2  IS  A  LINEAR  FUNCTION  OF  L  WITH  A  SLOPE  OF 
4ir2C.  (b)  THE  LINEAR  RELATION  BETWEEN  T2  AND  L  FOR  A 
1/8  INCH  POLYSTYRENE  SPECIMEN. 
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FIGURE  6.  CURVES  ILLUSTRATING  THE  BUILD  UP  OF  ELECTRIC  CHARGE  ON  THE  PLATES 
SURROUNDING  THE  POLYSTYRENE  SPECIMENS.  VALUES  FOR  THE  CHARGE,  Q,  WERE 
OBTAINED  FROM  THE  TIME  INTEGRAL  OF  THE  CURRENT  THROUGH  THE  LOAD  RESISTOR. 
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